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Dissipation effects in percolating quantum Ising magnets 
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(— i Abstract 
I 

(ZJ We study the effects of dissipation on a randomly dilute transverse-field Ising magnet at and close to the percolation threshold. 
'r^^ 5 , ° r weak transverse fields, a novel percolation quantum phase transition separates a super-paramagnetic cluster phase from an 
^ ♦ inhomogeneously ordered ferromagnetic phase. The properties of this transition are dominated by large frozen and slowly fluctuating 
1^3 percolation clusters. Implementing numerically a strong-disorder real space renormalization group technique, we compute the 
C"j low-energy density of states which is found to be in good agreement with the analytical prediction. 
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In dilute quantum magnets, the combination of geomet- 
ric and quantum fluctuations can lead to unconventional 
low-temperature properties such as singular thermody- 
namic quantities in quantum Griffiths phases as well as 
exotic scaling at the quantum phase transitions pQ. In 
many real systems, the magnetic degrees of freedom are 
coupled to an environment of "heat bath" modes. The 
dissipation caused by the bath is known to qualitatively 
change the properties even of a single quantum spin [2]. 
Here, we show that even small dissipation dramatically 
changes phases and phase transitions of a randomly dilute 
quantum Ising magnet. 

We consider the rf-dimensional site-dilute random 
iransverse-field Ising model given by the Hamiltonian 



> Hi = - 

X- 
S3 



(i,j) 



J, 



(i) 



where the Pauli matrices erf and of represent the spin com- 
ponents at site i, the exchange interaction JJy > couples 
nearest neighbor sites, and the transverse field hi controls 
the quantum fluctuations. Dilution is introduced via ran- 
dom variables K t which can take the values and 1 with 
probabilities p and 1 — p, respectively. We now couple each 
spin to a local ohmic bath of harmonic oscillators, 



+ 2 X i^ a i( a i,n + a i,n) 



(2) 
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i,.„ imii a\ are the annihilation and creation oper- 
ators of the n-th oscillator coupled to spin i; v^ n is its fre- 
quency, and \i yTl is the coupling constant. The baths have 



spectral functions £i(u>) = TEn^in^( w 



Vi 



2naiUje~ u '/ ni with on the dimensionless dissipation 
strength and tti the cutoff energy. 

For small transverse fields in the absence of dissipation 
(pti — 0), this system undergoes an unconventional phase 
transition controlled by an infinite-randomness fixed-point 
(p = Pc) from a ferromagnetic phase (p < p c ) to a quantum 
Griffiths paramagnetic phase (p > p c ), where p c is the ge- 
ometric percolation threshold. Close to the transition [5], 
the low-energy density of states (DOS) p (e) is dominated 
by large clusters tunnelling with an exponentially small 
frequency, A s w /iexp{— Bs}, where s is the cluster size, 
B = In (J/h), and J and h are the typical values of Jij and 
hi, respectively. Averaging over the cluster size distribution 
n s [1] yields (in the paramagnetic phase) 



' n s 6 (e - A s ) 



1 



(3) 



where z ~ £ D is the dynamical exponent, £ is the geomet- 
ric correlation length and D is the fractal dimension of the 
infinite cluster. n s (£) = s~ Tc f (ts ac ), where t = p — p c , r c 
and cr c are classical percolation exponents and / is a uni- 
versal function: f(x) ~ { — (cix) 1 /°' c }, for x > Oand/(x) ~ 
exp{— (c2xy i ~ 1 / d ^ rTc }, for x < 0, with c\ and C2 being 
nonuniversal constants of order unity [4]. At the percola- 
tion threshold, z is formally infinity. In this case, because 
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geometrical fluctuations are maximal, 
1 



P(e) 



eln Tc (/i/e) 



(4) 



becomes even more singular. Below p c , the infinite cluster 
contributes with a delta peak at zero energy to the DOS. 

Dissipation dramatically changes the above scenario (for 
instance, consider the case where on = a and Qi = Q). 
It was shown that for small fields the effective dissipa- 
tion strength is additive [5], i.e., the effective dissipation 
strength of a cluster of size s is a s = as. Therefore, a clus- 
ter of size s tunnels with a renormalized frequency A s ~ 
A.CA./n) "/^- *) = hexp{-bs/(l - as}}, if a, < 1. If 
a s > 1, the cluster becomes static, i.e., A s = [2]. The 
constant b = B + aln(Tl/h). Note that dissipation strongly 
suppresses the tunnelling of large clusters of sizes near s c = 
1/a. This qualitatively changes the physical behavior in 
comparison to the dissipationless case. The paramagnetic 
quantum Griffiths phase is replaced at low-energies by a 
classical super-paramagnet phase in which large clusters 
behave nearly classically [5] . 

Computing the low-energy DOS, we find that 



S<S C 



1 



eh/ (h/e)' 



(5) 



with 4> = 2, for energies e < e cross , where e croS s = 
hexp {— b/a} is a crossover energy scale above which weak 
dissipation becomes unimportant and p (e) reduces to 
Eq. ([3]) or (|4j) . Furthermore, the singular dependence of 
the DOS Eq. (J5]) does not depend on the p, meaning it is 
the same below, at and above to the percolation threshold. 
The infinite cluster (when p < p c ) and the frozen clusters 
(when s > s c ) contribute to a delta peak to the DOS. 
Thermodynamical quantities in the above scenario show 
unconventional behavior as discussed in Ref. 5J . 

We now use a generalization of the strong-disorder 
renormalization group method in order to compute the 
DOS [617] . We search for the strongest bond in the system 
max{/ii, Jij}. (a) In the case of a coupling, say J23, spins a 2 
and 03 become locked in a ferromagnetic state, and thus, 
are replaced by a single effective spin a 2 ■ The renormalized 
transverse field reads h 2 — h 2 h 3 (J 2 /Q 2 )° 2 (^/f^)" 3 / Ji- 
The effective dissipation strength is a 2 = a 2 + a 3 and the 
new cutoff Cl 2 = J 2 . The coupling between the new spin 
and spin cr* is J 2 i = J 2 i + J^i- (b) If the strongest bond is a 
field, say h 2 , then we proceed as follows, (b.i) If Vl 2 > h 2 , 
we integrate out the oscillators until Q 2 = h 2 . Hence, 
we renormalize h 2 — h 2 (h 2 /Q 2 ) a2 and Q 2 — h 2 . (b.ii) If 
fl 2 = h 2 , we freeze the spin (T2 in the direction of the trans- 
verse field h 2 . Then, we remove it from the system and 
renormalize the remaining couplings = Jy + J 2 iJ 2 j/h 2 . 

Iterating these steps, the DOS is obtained by building the 
histogram of effective tunneling rates e of the percolation 
clusters. Fig. [1] shows the distribution P(x), where x 
r: 1 < . . for p — p c and p — 0.5 for a two dimensional square 
lattice of size LxL, and coupling constants Jij (transverse 



fields hi) uniformly distributed between and 1 (0 and 0.1). 
According to Eq. ([5]), P(x) ~ x~^, with <f> — 2 regardless 
the value of p. As shown by the fits in Fig. Q] the data 
are in very good agreement with this analytical prediction. 
Moreover, for p = 0.5 7^ p c , the data disagree with the 
prediction of the undamped system Eq. ([3]), showing the 
importance of the dissipation for the low-energy DOS. 
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1. Distribution of low-energy tunneling rates (a) for 
p c = 0.407253 and (b) p = 0.5 for different system sizes L. The 



distribution were built decimating 2 000 systems for L 
100, and 8 000 systems for L = 50. 
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In conclusion, we have shown that dilute quantum Ising 
magnets with Ohmic dissipation undergo a novel perco- 
lation quantum phase transition between a ferromagnetic 
and a classical super-paramagnetic phase. The low-energy 
DOS displays striking singularities due to the suppression 
of tunneling by damping. Using a strong-disorder renormal- 
ization group method, we have computed the low-energy 
DOS numerically and found it to be in good agreement 
with the analytical prediction. 
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